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ABSTRACT 


A^aodel  of  optical  and  infrared  wave  propagation  in  scattering 


media  is  derived  from  Maxwell's  equations.  It  allows  for  the  calcula- 


,v 

ft 


tlon  of  the  irradiance  profiles  of  the  forward-  and  the  backward- 
propagating  beams  which  result  from  the  interaction  of  a  coherent  beam 
with  a  scattering  medium.  The  method  accounts  for  all  orders  of  scat¬ 
tering  and  is  applicable  to  inhomogeneous  media.  The  range  of  validity 
is  limited  by  the  condition  that  the  fractional  volume  occupied  by  the 
particles  to  power  must  be  much  smaller  than  unity.  {■''.*  ,  n  rioj, 

c  i/*rtk 


A 


RESUME 


On  a  mis  au  point,  i  partir  des  Aquations  de  Maxwell,  un  module 
de  propagation  des  ondes  optiques  et  infrarouges  dans  un  milieu  dlffu- 
sant.  Ce  module  permet  de  calculer  les  prof 11s  de  1' Intensity  lumi- 
neuse  des  faisceaux  se  propageant  vers  l'avant  et  vers  l'arriire  et  qui 
rAsultent  de  1' interaction  d'un  falsceau  coherent  avec  les  partlcules. 
La  mCthode  tlent  compte  de  tous  les  ordres  de  diffusion  et  peut  s'ap- 
pllquer  aux  milieux  InhomogAnes.  Le  domalne  de  valldltA  est  limit A  par 
la  condition  sulvante:  la  fraction  du  volume  occupA  par  les  partlcules 
i  la  puissance  i  dolt  fttre  beaucoup  plus  petite  que  l'unlti. 
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1.0  INTRODUCTION 


A  characteristic  feature  of  modern  warfare  Is  the  widespread  use 
of  electro-optical  systems  designed  to  enhance  the  effectiveness  of 
various  weapons.  However,  poor  weather  conditions  and/or  artificial 
obscurants  can  severely  degrade  their  performance.  Hence,  a  require¬ 
ment  exists  for  understanding  and  predicting  the  propagation  of  elec¬ 
tromagnetic  waves  under  adverse  conditions  created  mostly  by  the  pres¬ 
ence  of  suspended  natural  or  man-made  aerosols. 

Propagation  in  thin  or  tenuous  clouds  is  well  understood  and 
documented  (Refs.  1  and  2).  It  is  governed  by  simple  and  exact  single- 
scattering  equations  which  explain  numerous  atmospheric  phenomena. 
However,  for  military  applications,  there  are  many  instances  of  low 
visibility  conditions  where  the  single  scattering  hypothesis  fails.  In 
such  dense  clouds,  the  optical  rays  undergo  many  scattering  events 
before  escaping  the  medium,  reaching  a  target,  or  being  detected. 

These  multiple  scatterings  have  nonnegliglble  and  often  dominant 
effects  on  the  propagated  wave.  Analytic  methods  (e.g.  Refs.  3-5), 
transport  methods  (e.g.  Refs.  6-15),  and  Honte  Carlo  simulations  (e.g. 
Refs.  16-17)  have  all  been  used  to  study  this  problem.  The  mathematics 
are  complex  and  no  general  solution  has  been  obtained  yet.  The  numeri¬ 
cal  analysis  based  on  the  Monte  Carlo  method  is  time-consuming  and 
provides  empirical  relations  only.  Fortunately,  there  are  special  and 
practical  cases  where  simplifications  are  possible,  which  result  in 
useful  approximate  theories.  The  most  accurate  of  these  approaches  was 
developed  at  Defence  Research  Establishment  Valcartier  (DREV)  (Refs. 
11-15).  It  consists  of  the  exact  series  solution  of  the  small-angle 
approximation  to  the  radiative  transfer  equation.  The  model  is  very 
convenient  for  numerical  computations  and  gives  valuable  insight  into 
the  relative  effects  of  the  various  scattering  orders.  However,  it  is 
restricted  to  homogeneous  media  and  it  does  not  appear  readily  applica¬ 
ble  to  the  important  inverse  problem  of  interpreting  lidar  measurements 
for  the  determination  of  the  optical  properties  of  the  medium. 
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A  model  based  on  a  different  theoretical  approach  is  developed 
here*  It  steins  from  the  Maxwell's  equations  of  electromagnetic  wave 
propagation  and  not  from  the  radiative  transfer  equation.  The  method 
predicts  the  irradlance  profiles  of  the  forward-  and  backscattered 
beams  including  the  Ildar  configuration.  All  orders  of  scattering  are 
accounted  for  but  they  are  not  individually  identified.  The  theoreti¬ 
cal  basis  is  outlined  in  Chapter  2.0.  Chapter  3.0  describes  the  sto¬ 
chastic  approach.  The  original  features  of  the  model  are  introduced 
and  discussed  in  Chapter  4.0,  and  Chapter  5.0  gives  the  resulting  prop¬ 
agation  equations.  Sample  solutions  for  beam  waves  and  plane  waves  are 
derived  and  analysed  in  Chapters  6.0  and  7.0  respectively. 

This  work  was  performed  at  DREV  between  January  and  December 
1983  under  PCN  33B07,  Atmospheric  Propagation  of  Laser  Beams* 

2.0  THEORETICAL  BACKGROUND 


The  object  of  this  work  is  to  describe  the  propagation  of  the 
forward-  and  backward-going  electromagnetic  waves  in  a  particulate 
medium.  We  restrict  the  analysis  to  situations  where  polarization 
effects  are  negligible.  Under  this  approximation,  the  wave  electric 
field  E  can  be  considered  scalar  and,  from  Maxwell's  equations,  it 
satisfies  the  equation 

[i] 

at2  c2 

where  V2  is  the  three-dimensional  Laplacian  operator,  t  is  the  time,  p' 
is  the  relative  magnetic  permeability,  e'  is  the  relative  dielectric 
constant  of  the  medium,  and  c  is  the  speed  of  light  in  free  space.  The 
right-hand-side  function  S  is  the  source-sink  term  arising  from  the 
Interaction  of  the  forward-  and  backward-propagating  waves  through  the 
particles  in  the  medium.  S  is  related  to  the  current  and  charge  densi¬ 
ties  induced  on  or  within  the  particles.  However,  it  will  not  be 
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formulated  In  this  exact  fashion  here;  It  will  be  modeled  as  described 
later  In  Section  4.1. 

We  consider  only  one  temporal  Fourier  component  of  the  electric 
field.  If  u  denotes  the  angular  frequency  of  the  wave,  the  temporal 
Fourier  transform  of  eq.  1  yields 

V2g  .  Ifelella;2  E  .  s.  [2] 

c2 

where  the  tilde  (~)  denotes  a  Fourier  transform.  Equation  2  was  sim¬ 
plified  through  the  approximation  w  >>  (p's')"1  which  Is  well 

justified  for  optical  and  infrared  waves.  The  product  (p' e')  Is  the 
square  of  the  complex  refractive  Index  n  and  it  Is  written 

p'  e'  -  n2.  [3] 

We  seek  a  solution  of  the  form 

E  -  A  exp  [lk(  ±*  +  ♦)  ],  [4  ] 

where  z  Is  the  coordinate  along  the  main  direction  of  propagation, 

1  ■  k2  *  u2(n  2  -  n,  2)/c2  is  the  wave  number,  n  and  n.  are 

ro  lo  ro  lo 

respectively  the  real  and  imaginary  parts  of  the  reference  refractive 
Index,  and  A  and  $  are  functions  to  be  specified  below.  On  substltut- 
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aA, 


dz  +  .  XAf  +  *  Af  Z  .  Lf  '  ~  “f  2k  *  "f 


V,  +  kW  A,  -  -ir-  72A. 


Sf  e‘ik(z+^f>. 


[7] 


z\  +  *\z 


Zb 


+  kWA. 

D 


Sb  e~ik(z+*b\[8] 


where  the  subscripts  f  and  b  refer  to  the  forward  and  the  backward  wave 
respectively,  and  where  E  7$f  and  E  74^.  The  functions  U  and  W  are 
defined  as  follows: 


U 


ro 


+  n 


lo 


[9] 


W  - 


nrni 


[10] 


where  nf  and  n^  are  the  real  and  Imaginary  parts  of  the  instantaneous 
and  local  complex  refractive  Index  of  the  particulate  medium. 


Two  types  of  separation  were  performed  In  obtaining  eqs.  5-8; 
first,  the  separation  Into  the  forward  and  the  backward  waves  which  has 
its  basis  in  the  +  and  -  signs  in  eq.  4,  and  second,  the  separation 
into  an  equation  for  the  phase  front  angle  V  and  the  amplitude  A.  The 
latter  separation  is  different  from  the  more  common  practice  of  factor¬ 
ing  out  the  real  and  imaginary  parts  in  eq.  2.  The  basis  is  to  uncou¬ 
ple  the  equation  of  V  (or  $)  from  that  of  the  complex  amplitude  A. 

Equations  5  and  6  are  the  eikonal  equations,  or  the  geometric 
optics  equations,  of  the  forward  and  backward  waves.  The  surfaces 
(z+$)  ■  constants  are  the  geometric  phase  fronts  and  the  vector  (V  +  k) 
gives  the  direction  of  the  optical  rays,  where  k  is  the  unit  vector 
along  the  z-axis.  Hence,  from  eqs.  4,  7  and  8,  it  follows  that  A  is 
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the  complex  amplitude  defined  on  the  geometric  phase  front.  This  com¬ 
plex  amplitude  embodies  the  phase  perturbations  induced  by  diffraction. 

3 . 0  STOCHASTIC  APPROACH 


For  the  present  analysis,  we  model  the  particulate  medium  as  a 
space-time  random  function  of  the  complex  index  of  refraction  n.  This 
function  can  be  specified  from  the  given  properties  of  the  medium,  l.e. 
the  number  density  of  the  particles,  their  shape,  their  size  distribu¬ 
tion,  the  distributed  complex  refractive  index  within  the  particles, 
and  the  complex  refractive  index  of  the  surrounding  medium.  Hence,  the 
quantities  V  and  A  are  random  functions  and  the  governing  eqs.  5-8  are 
stochastic  equations. 

There  is  no  known  general  method  of  solving  for  the  random  func¬ 
tions  V  and  A.  In  any  case,  this  would  yield  much  more  Information 
than  is  required  in  practice.  The  quantities  of  interest  are  the  aver¬ 
ages  and  one  possible  approach  consists  in  using  the  stochastic  eqs. 

5-8  to  derive  deterministic  equations  for  these  averages  or  statistical 
moments.  In  the  present  application,  we  are  interested  in  the  moments 
<A>  and  <AA*>  which  are  respectively  the  average  field  amplitude  and 
the  average  irradlance. 

The  method  of  derivation  is  as  follows.  First,  the  functions  U 
and  U  representing  the  medium  properties  are  written  as  sums  of  an 
average  and  a  random  part,  l.e. 

U  -  <U>  +  u  ;  <u>  -  0,  [ll] 

W  "  <W>  +  w  ;  <w>  “  0,  [l2] 


where  the  pointed  brackets  denote  ensemble  averaging.  The  same  decom¬ 
position  is  performed  for  the  dependent  wave  front  function 
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V  -  <V>  +  v  ;  <v>  -  0.  [13] 

The  field  amplitude  A  Is  separated  In  a  slightly  different  way  Into  a 
coherent  (subscript  c)  and  an  incoherent  (lower  case  letter)  contribu¬ 
tion 


A  -  Afi  +  a  ;  <a>  ■  0.  [14] 

The  component  a  is  random  and  represents  the  scattered  amplitude. 

The  component  A is  not  fully  deterministic.  There  are  randomly 

distributed  holes  in  the  profile  of  A  left  by  the  presence  of  the 

c 

scatterers. 


If  we  substitute  eqs.  12-14  for  W,  V  and  A^  in  eq.  7  and  drop 
the  subscript  f  for  brevity,  we  obtain 


+  £  +  <V>  .  7A  +  <V>  .  7a  +  v  .  7A  +  v  .  7a 

oz  oz  ~  ~  c  ~  ~  ~  ~  c  ~  ~ 

+  i  A  V  .  <V>  +  -i-  a  7  .  <V>  +  4  A  7.v+— a7.v 

2  c  ~  2  2  c  ~  ^  2  <v  m 


+  k  <W>  A  +  k  <W>  a  +  kw  A  +kwa 
c  c 


-  L  y2  A  _  y2  -  3f  e”ik^^ 

2k  V  c  2k  V  *  S  * 


[15] 


The  terms  of  the  form  V  .  VA  are  transport  terms.  They  account  for  the 
propagation  of  the  amplitude  A  along  the  rays  of  direction  (krf-V) .  The 
terms  y  A  7  .  7  represent  the  compressibility  of  the  amplitude.  A 
increases  if  the  rays  converge  and  decreases  if  they  diverge.  Accord¬ 
ing  to  this  general  interpretation,  we  can  set  in  eq.  15 
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which  states  that  the  coherent  amplitude  Is  not  propagated  by  the  inco¬ 
herent  or  scattered  rays.  Physically,  this  is  explained  by  the  fact 
that  »  0  immediately  past  the  particles  so  that  it  cannot  be  trans¬ 
ported  by  the  scattered  rays  that  originate  from  the  particles.  On  the 
other  hand,  the  term  A£  7  .  v  is  not  zero  instantaneously  since  it 
models  the  action  of  the  scattering  rays  on  the  nonzero  coherent  ampli¬ 
tude  that  impinges  on  the  particles.  The  term  y  7  .  v  is  responsi¬ 
ble  for  the  depletion  of  the  coherent  amplitude,  l.e.  its  transforma¬ 
tion  into  scattered  amplitude.  It  can  be  rewritten 

4  A  V.v«i<A>V»v+i(A-^l>)  V.v  *  [l7  ] 

Taking  the  ensemble  average  of  eq.  17,  we  have 

4  <A  V  .  v>  -  4  <(A  -  <A  >)  V  .  v  >.  [l8  ] 

2  c  2  c  c  '  ~  lj 

Whereas  the  function  (A  -<A  >)  is  characterized  by  narrow  negative 

c  c 

spikes  created  by  the  shadows  of  the  particles  and  does  not  depend  at 
all  on  their  refractive  index,  the  random  phase-front  function  V.v  is 
distributed  more  uniformly  over  the  complete  space  and  is  a  strong 
function  of  the  refractive  index  of  the  particles.  Therefore, 

(Ac~<Ac>)  is,  for  all  practical  purposes,  not  correlated  to  V.v  and  it 
follows  from  eq.  18  that 

J  «£  I  '  I’  *  0  .  [»] 

and  similarly  that 

<Afi  w>  *  0.  [20] 

Using  eqs.  9-20  and  applying  them  to  both  the  forward  and  the 
backward  fields,  noting  that  <A>*^Ac>,  dropping  the  subscript  c  when 
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not  necessary,  and  taking  the  ensemble  average  (term  by  term)  of  the 
stochastic  eqs.  5-8,  we  derive  the  following  set  of  deterministic  equa¬ 
tions  for  the  average  wave  front  angles  and  wave  amplitudes: 


8<Ff> 


s r  +  y  •  l  y  •  l  [<<'>  -  •  vl* 


»y 


+  y  •  l  y '  Z  t<0>  -  •  v)« 


9^Af  > 


+  <V>  .  V  <Af>  +  |  <Af>  V  .  <Vf>  +  k  <W>  <Af 


“  k  v2<V  “  “k  <wa£>  ~  l  •  W  +  *  ^fZ  •  V 

"lk<*  +  *f> 

+  <Sf«  >. 

a  <*.  > 

-  -sr  +  y  •  z  y +  *  v  s  •  +  k  y 

"  Ik  v2<Ab>"  “k  <WV  ■  X  •  Vb' +  *  •  V 


+  <Sbe’lk(’*+V>. 


[21] 

[22] 


[23] 


[24] 


Subtracting  eqs*  23  and  24  from  eqs*  7  and  8  respectively,  and  using 
eqs*  9-20,  we  arrive  at  the  equations  for  the  random  amplitudes  a^  and 

V 
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"  .  Vaf+|af  £  •  -  £  .  %  «f>  +  y  <af  £  .  v^} 


+  [Sf  exp[-ik(z+$f)  ]  -  <Sf  exp[-lk(z+$f)]>]. 


‘?+V*X8b  +  iabi‘ V  +  k<w>ab'k72ib 


-  -  {kwas  -  k<waK>}  -  kwAcb-  |  A cb  •  £  ^ 


"  feb  *  l  ab  +  1  ab  l  *  l  *  %  V :  +  1  <ab  l  •  V* 


+  {Sb  exp[-ik(-z+^b)  ]  -  <Sfe  exp[-ik(-z+4>b)]>}.  [26] 

Finally,  the  equations  for  the  average  forvard-  and  backward-scattered 
lrradlance  are  readily  obtained  from  eqs.  25-26  and  16-20.  They  are 


a  > 


+  <V_>  .  V  <tt.*5  >  +  <aea*  >  V  .  <V,>  +  2k  <W>  <a_«*  > 
"t  ~  f  f  f  f  ~  ff 


-  ik  X  •  [<■?  l  v  -  l  af>] 


-  2k  <waja*>  -  k  <wa^>  <A^>*  -  k  <waj>  <A^> 

-  i  <A#>  <a*  7  .  v  >  -  ~  <A  >*  <a,  7  .  v  > 

2  f  r  ~  2  f  f  ~  ^ 

-  v  •  +  <*£  Sf  exp  [-lk(z+7f )  ]> 

+  <af  S*  exp[+lk(z+*f)]>, 


[27] 
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-k  l  [<*{  IV  VH’l 

■  -  2k  <wa^a*>  -  k  <Va^>  <A^>*  -  k  <wa*>  <A^> 

-  j  <Ab>  <a*  V  .  vfc>  -  y  <Ah>*  <afe  7  .  vfc> 

-  £  •  %abab>  +  <•£  Sb  exp[-ik(-z+^)]> 

+  ^  S*  expf+lkC-z+^h)  j>.  [28] 

The  system  of  deterministic  eqs.  21-24  and  27-28  is  mathemati¬ 
cally  unclosed  insofar  as  it  contains  more  unknowns  than  equations , 

such  as  the  moments  <v.a_>,  <a-7.v_>,  etc*  This  is  the  classic  closure 

if  f  ~ 

problem  always  encountered  in  the  treatment  of  statistical  phenomena 
governed  by  nonlinear  and/or  quasi-llnear  stochastic  equations,  such  as 
eqs.  5-8.  The  complete  mathematical  model  contains  an  infinite  set  of 
equations.  Hence,  workable  models  require  closing  this  hierarchy  of 
equations  at  a  practical  order. 

The  set  of  eqs.  21-24  and  27-28  is  exact  in  the  framework  of 
unpolarized  forward-  and  backward-propagating  waves  scattered  by  a 
particulate  medium,  but  it  is  indeterminate.  The  essence  of  the  pro¬ 
posed  model  will  consist  in  deriving  constitutive  relations  to  relate 
the  unknown  statistical  moments  on  the  right-hand  sides  of  eqs.  21-24 
and  27-28  to  the  lower  order  functions  <A>  and  <aa*>. 
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4.0  MODEL 

4.1  Source-Sink  Terms 

The  first  element  of  the  model  Is  the  formulation  of  the  source- 
sink  terms.  These  terms  are  responsible  for  the  interchange  of  energy 
between  the  forward-  and  the  backward-propagating  waves.  This  interac¬ 
tion  arises  because  of  the  backward  reflection  and  scattering  by  the 
individual  particles  of  a  fraction  of  the  impinging  radiation.  Instead 
of  trying  to  solve  this  problem  by  rigourously  applying  the  electromag¬ 
netic  boundary  conditions  at  the  interface  of  the  particles  and  the 
surrounding  medium,  we  model  the  backscattering  source-sink  terms  in  a 
form  mathematically  similar  to  the  forward  scattering  interaction  as 
expressed  by  the  terms  of  the  form  -j  AV.V  of  eq.  IS.  Hence,  we  write 

’  7  ‘lk<*+^)  I  •  -  Zt' >•  [2»] 

h  -  T  .lk<'**V  Af  7  .  (V-  -  J^),  [30] 

where  v^  and  are  the  backscattered  random  geometric  phase  front 
angles  of  the  forward  and  backward  wave  respectively.  They  will  be 
determined  by  application  of  the  reflection-refraction  laws  at  the 
interface  of  the  particles. 

The  mechanism  of  separation  into  a  forward  and  a  backward  compo¬ 
nent  modeled  by  eqs.  29  and  30  is  somewhat  arbitrary.  The  original 
Maxwell's  equations  contain  no  such  distinction.  In  effect,  eqs.  29 
and  30  amount  to  a  special  and  simplified  treatment  of  the  boundary 
conditions  at  the  interface  of  the  particles.  Although  heuristic,  this 
separation  into  forward-  and  backscattering  is  very  convenient  in  prac¬ 
tical  applications. 
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4.2  Hypotheses 

The  second  most  important  element  of  the  model  Is  the  derivation 
of  the  constitutive  relations  for  the  higher  order  unknown  moments. 

The  derivation  is  based  on  the  following  hypotheses: 

Hypothesis  1.  The  random  phase  front  angle  v  and  amplitude  a  are 
only  weakly  correlated.  This  approximation  is  consistent  with  the 
governing  eqs.  5-8,  for  they  imply  that  v  is  independent  of  a  and  that 
the  random  a  is  the  result  of  repeated  interactions  with  the  phase 
front  angle  v  over  the  complete  propagation  path  and  not  only  at  the 
point  of  observation.  What  the  hypothesis  means,  for  instance,  is 
that 

<va>  «  <v  v  <aa*>^. 


[31] 


Such  a  relation  was  experimentally  verified  in  the  random  continuous 
medium  of  turbulence  (Ref.  18).  In  particulate  media,  we  have  no  such 
direct  measurements  but  the  hypothesis  is  self-consistent.  For  exam¬ 
ple,  the  result  obtained  for  <va >  satisfies  the  Inequality  [31  ]  if 
<v.v>^  «  1.  From  the  equation  for  v,  we  find  that  the  standard  devia- 
tlon  of  v,  i.e.  <v.v>  ,  is  of  the  order  of  the  fraction  of  volume 
occupied  by  the  particles  to  the  power  i.  For  naturally  occurring 
aerosols,  this  quantity  rarely  exceeds  10~5  which  gives  <v.v>^  »  0.06. 
Hence,  the  results  derived  from  Hypothesis  1  are  self-consistent  under 
most  practical  conditions. 


Hypothesis  2.  The  cross  correlations  between  quantities  per¬ 
taining  to  the  forward  and  backward  waves  are  negligible,  e.g. 
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Although  there  Is  some  interdependence  between  the  forward  and  backward 
waves,  the  local  quantities  depend  mostly  on  the  random  properties  of 
statistically  independent  regions  of  the  scattering  medium.  Hence,  the 
hypothesis  appears  well  justified,  especially  since  the  cross-correla¬ 
tion  terms  are  always  neglected  by  comparison  with  direct-correlation 
term s. 

Hypothesis  3.  Although  not  essential,  we  make  the  paraxial 
approximation.  This  considerably  simplifies  the  mathematics  and  fol¬ 
lows  from  our  finding  that  the  standard  deviation  of  the  random  phase 
front  angle  is  not  larger  than  <'0.06  for  most  practical  applications. 

Hypothesis  4.  The  particulate  medium  is  statistically  homoge¬ 
neous  and  isotropic  over  a  domain  the  size  of  the  beam  diameter.  This 
approximation  is  well  justified  for  propagation  of  low-divergence  laser 
beams. 


Hypothesis  5.  The  random  phase  front  angle  is  statistically 
homogeneous  and  isotropic  in  the  plane  transverse  to  the  main  direction 
of  propagation,  the  0-z  axis.  This  is  consistent  with  Hypotheses  3  and 
4  which  guarantee  that  all  the  rays  reaching  a  plane  z  traverse  statis¬ 
tically  equivalent  paths. 

Hypothesis  6.  The  complex  amplitude  covariance  function  is 
quasl-homogeneous  and  quasl-lsotroplc  in  the  transverse  plane.  More 
specifically,  we  assume  that  it  is  of  the  form 

<a^*l»£l)a*(z2»l2)>  ■  F[zl»*2J<*l’**2>^2]  *  f  tz  l,z 2'  S.rS.2]*  t33  ] 

where  F  and  f  are  unspecified  functions  used  to  illustrate  the 
functional  dependence  only.  This  is  a  standard  approximation  made  in 
the  presence  of  an  average  gradient.  It  is  not  essential  but  very 
convenient. 
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4.3  Constitutive  Relations 


From  Hypotheses  3-5,  eqs*  21  and  22  simplify  to 


a<vf> 

aS 


+  <v  > 

~f 


o. 


[34] 


a<v  > 

~o 
dz 

where  denotes  the  two-dimensional  gradient  operator  In  the  trans¬ 
verse  plane,  l.e.  with  respect  to  the  two-dimensional  position  vector 
r .  Equations  34  and  35  show  that  the  average  phase  front  angle  is 
independent  of  the  statistical  properties  of  the  medium.  The  solution 
of  eq.  34  for  an  initially  spherical  phase  front  with  radius  of  curva¬ 
ture  (focal  distance)  F  is  given  by 


+  <V  > 
~b 


V.<V> 


0, 


[35] 


<!*>  -  [36] 

From  eq.  27  and  the  assumption  of  statistical  symmetry  of  the  scat¬ 
tering  pattern  with  respect  to  the  incident  ray  it  follows  that 

•V  '  -%*  '  -£/<*-»>•  [33] 

The  first  step  toward  the  determination  of  the  constitutive 
relations  is  the  solution  for  the  random  amplitudes  a^  and  a^.  Solving 
implicitly  eqs.  25  and  26  with  the  respective  boundary  conditions 
af(z-0)  -  0  and  afe(z«Z)  -  0,  we  find 

2  2  oo 

af(z'0  exp  W*^***'}  I  ^  II  d2» 

o  o  — • 


exp  {• 


1k 

2(z-x) 


F-z  «  F-x 
F-x  ‘  F-z 


r  -  s 


o 


k<W>dz  * } ,  [38] 
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Z  Z  00 

ab(z»~)  “  TJ  exp  W  k<W>dz'}  /  //  d28  gb(x,8) 

z  z  -» 

•  exp  ^2(^0  f=£  I  E  -  £  j2  +  /  k<W>dz '  } ,  [39] 

X 

where  g,(x,s)  and  g  (x,s)  stand  for  the  right-hand  sides  of  eqs.  25  and 
26  respectively  and  where  the  solutions  [36]  and  [37]  have  been  substi¬ 
tuted  for  <V^>  and  <]£b>*  Equations  38  and  39  are  Implicit  solutions 
for  the  forward  and  backward  random  amplitudes  since  the  functions  g^ 

and  g,  contain  a,  and  a,  . 
d  r  b 

The  unknown  higher  order  moments  are  <wa,>,  <v.a,>,  <a,V,  •  v_>, 

r  ~r  r  f~L 

etc*  The  procedure  to  relate  these  moments  to  the  lower  order  <A>  and/ 
or  <aa*>  is  basically  the  same  for  all  of  them.  The  implicit  solution 
for  a^  (or  a^)  is  first  multiplied  by  w,  or  v^,  or  V  .  v^,  etc...,  and 
the  products  are  then  ensemble  averaged.  This  produces  many  terms 
under  the  Integral  operator  of  eq.  38  or  39.  Most  of  them  can  be  ne¬ 
glected  in  relation  to  a  leading  term  on  account  of  Hypothesis  1,  and 
further  mathematical  simplifications  can  be  worked  out  using  Hypotheses 
2-6.  These  derivation  steps  are  lengthy  and  will  not  be  repeated  here; 


however  they  are  similar  to  those  followed  in  Refs.  19  and  20. 

resulting  constitutive  closure  relations  are 

The 

k 

<wa^>  a 

l 

‘  2  aaf 

<Af>, 

[40] 

k 

<wa.  >  * 

D 

1 

2  °ab 

t*l] 

k 

<wa^a^> 

1  _ 

*  -  |  Re 

^aafs  • 

]  <afa*>» 

[42] 

k 

<wa, a*> 

D  D 

1  _ 

*  -  j  Re 

t°ab8  ■ 

)  <vS>- 

[43] 

'  '-'.v  >'  -v  .v 


•LVO  vW.rv'TV 
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Zl  '  Xf>  "  '  °cf  >* 

h  *  V  *  "  °cb  ^b'* 

ik<VV 


ik(*f-<^) 


ik(^-*f) 


ab  Z±  •  <s  -  -  ‘Gf<V* 


af  Xi  •  ■  Xb>>  3  -  ^b'V* 


a*  ab  Xi  •  (tb“  £>>  +  c-c- 


*  2  Ke[a;b]  <ab  afe*>  -  2  Re[tTf  ]  <af  af*>. 


Sf J  f  “f 


%  ®f>  *  °» 


%  V  3  °* 


%  afaf>  “  -  R®[Df  ]  •  XjL 

%  abab>  3  -  RetebJ  •  ii 


[44] 

[45] 

[46] 

[47] 


[48] 

[49] 

[50] 

[51] 

[52] 


where  c.c.  stands  for  the  complex  conjugate  of  the  preceding  term  and 
where  the  quantities  underlined  by  a  double  tilde  are  two-dimensional 
tensors  or  dyadics. 


There  are  several  constitutive  relations  but  they  can  all  be 
classified  in  one  of  three  categories:  the  particle  absorption  terms, 
eqs.  40-43;  the  aerosol  forward-  and  backscattering  terms,  eqs.  44-48; 
and  the  diffusion  terms,  eqs.  49-52. 

The  coefficients  of  eqs.  40-52  are  given  by  the  Integral  opera¬ 
tor  38  or  39  operating  on  various  covariance  functions.  To  simplify 
the  algebra,  we  restrict  applications  to  slightly  diverging  beams.  Ii 
that  case,  the  focal  length  F  of  eqs.  38  and  39  is  large  and  negative. 
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Moreover,  the  covariance  functions  operated  upon  have  longitudinal 
correlation  lengths  that  are  generally  much  smaller  than  | F— z |  and  thus 
the  major  contributions  to  the  integrals  come  from  values  of  x  in  the 
immediate  neighbourhood  of  z  compared  with  |F-z|.  Hence,  we  can  make 
the  approximation 


F-z 

F-x  “  1 » 


[53] 


and  the  following  simplified  operators  0+  and  0"  can  be  written  from 
eqs.  38  and  39: 

z  “"A  (z,  ® 

0+{h(z,r;Ti,£)}  -  -  ||  /  e  f  JJ  d2£  e  2n  h(z,r;n,£),  [54] 


Z-z  .  -A.(z,t)) 

0-{h(z,rjri,£)}  -  -  ^  ?  e 


2n  o  11 


•  lkp2 

//  d2£  e  2r)  h(z,r;ti,£),  [55] 


where  h(z,r;Ti, £)  is  a  dummy  function  and  where 


Af(z,Ti)  ■  /  k<W(z-x')>  dx', 

o 


A^( z , t))  -  /  k<W(z-x')>  dx'. 


[56] 

[57] 


In  terms  of  the  operators  0  and  0  ,  the  particle  absorption  coeffi¬ 
cients  are  given  by 


oaf  -  2k2  0+  {<w(l)w(2)>}. 

[58] 

aab  -  2k2  0-  [<w(l)w(2) >}, 

[59] 

aflfg  -  2k2  0+  (Gf(z,T),£)  <w(l)w(2)>}. 

[60] 

i 


M 
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aahs  "  2k2  °~  {Gh(*»Tl»£>  w(2>>}» 


the  scattering  coefficients  by 


ut  •  - 1  °f tXp  v  Vl>  V2)>>- 


°Cb  -  -  2  °' tep  v  V2> 


it  -  - 1  °~ i>t<2.''.6)  ip  v  [V2>  V2)> 


+  V‘1  £<*>>11. 


<b  -  -  7  °rf(W*.>be>  Vp  Vpl  [V«  £(2)> 


+  Vl>  Xi<2)>]}. 


®»f  *  '7  “ft'.l-E)  Xp  !p!  [%(!>  X^<2)> 


+  %  <»  i  <»>]!. 


°Ib  ■ "  i  °f{w*.'be)  Xp  V  i<2)> 
+  %  o>  i  <2>>ll. 


and  the  diffusion  coefficients  by 


8f  -  ?  0+  (Gf(2,%e)  <^(1)  X,  (2>>), 


M"'  ~f 


2b  -  i  °‘  K^e>  %(1>  Xb  [« 

where  (1)  and  (2)  refer  to  the  adjacent  points  (s,rf£/2)  and  (x-r|#r- 
£/2)  of  the  covariance  functions  and  where 


>  ■  - 


wammm 


vvv’.'VAW.'v'' 

rJ-  • .  \.\V- ,V*\ .Vis «SV-\>  >  YvVV t.' •  »\V -W.v\>WVvv*\  V .  vi^liUUiteVvUV. 


hHJWW 


UNCLASSIFIED 

19 


Gf(z,h.£)  ■  <«f(l)  «f*(2)>  /  <afaf*>. 

[70] 

-  <ab(D  «b*<2>>  /  <Vb*>. 

[71] 

H<z,rug)  -  <e  >. 

[72] 

The  list  of  33m bo Is  at  the  beginning  of  this  report  defines  the 
individual  absorption,  scattering  and  diffusion  coefficients. 


tatlon  Coefficients 


The  absorption  coefficients  0 ^  and  0^  depend  on  the  covariance 
function  4r(l)  w(2)>  which  involves  the  real  and  imaginary  parts  of  the 
random  refractive  index  as  defined  by  eqs.  10  and  12.  <w(l)  w(2)>  can 

be  determined  from  the  specified  medium  and  particle  properties,  l.e. 
the  complex  refractive  indices  of  the  medium  and  the  particles,  the 
shape  of  the  particles,  their  number  density,  and  their  size  distribu¬ 
tion.  Then,  oflf  and  0#b  can  be  calculated  by  direct  substitution  in 
eqs.  58  and  59.  Because  of  Hypothesis  4,  the  covariance  function  <w(l) 
w(2)>  is  independent  of  r  over  a  domain  the  size  of  the  beam  diameter. 
Hence,  the  coefficients  0flf  and  0^  are  functions  of  the  longitudinal 
coordinate  z  only.  These  coefficients  are  second-order  corrections  to 
the  absorption  coefficient  k<W>  of  eqs.  23-24  and  27-28. 


There  are  two  additional  absorption  coefficients  0^^  and 
which  depend  on  the  amplitude  correlation  functions  and  G^.  The 
latter  functions  are  not  known  but  their  correlation  lengths  are  ex¬ 
pected  to  be  much  greater  than  that  of  <w(l)  w(2)>  equal  to  the  average 

particle  size.  Hence,  both  G.  and  G.  may  be  replaced  (in  eqs.  60  and 

i  d 

67)  by  their  values  at  the  origin  which,  by  definition,  are  equal  to 
unity.  Therefore,  the  coefficients  oa^g  and  are  respectively 

equal  to  and  0^. 
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The  calculation  of  the  scattering  coefficients  (eqs.  62-67)  and 
the  diffusion  coefficients  (eqs.  68-69)  requires  knowledge  of  the  phase 
front  angle  covariance  functions.  These  cannot  be  inferred  directly 
from  the  medium  and  particle  properties  but  only  through  the  solution 
of  the  ray  propagation  equation.  The  procedure  is  to  derive  the  equa¬ 
tions  for  the  covariances  of  v,  and  v.  from  eqs.  5.  6,  21  and  22.  The 

~o 

result  for  slightly  diverging  beam  waves  and  approximately  Gaussian 

statistics  for  v,  and  v.  is 
~r  ~o 

%<D  vf(2)>  -  ^2  V1'  vf(2)>  -  i  V,  Vp^d)  ve(2)>: 

VD  tcm>>  -TSlplp  <a(1)  u(2)>1  [»] 

and  a  similar  equation  for  <v^(l)  v^(2)>.  The  source  term  <u(l)  u(2)> 
on  the  right-hand  side  of  eq.  73  can  be  determined  from  the  properties 
of  the  particulate  medium  in  the  same  fashion  as  for  the  function  <w(l) 
w(2)>  discussed  earlier.  Therefore,  the  solution  for  <v^(l)  v^(2)>  can 
in  principle  be  worked  out  from  eq.  73,  and  similarly  for  <vfe(l)  ^b(2)> 

The  solution  of  eq.  73  cannot  be  treated  exactly  because  of  the 
nonlinear  term.  No  satisfactory  method  of  solution  has  been  derived 
yet.  However,  a  preliminary  analysis  was  carried  out  and  order-of- 
magnltude  estimates  were  derived.  Since  the  analysis  is  Incomplete  and 
may  have  to  be  reexamined,  only  the  results  will  be  presented  here. 

First,  the  covariance  function  «v^(l)  y^(2)>  scales  as 

<v^.  Xf*  "  (7*1 

For  naturally  occurring  fogs,  we  obtain  from  eq.  74  that  the  standard 
deviation  of  the  forward  phase  front  angle  fluctuations  is  of  the  order 
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of  10  mrad,  In  qualitative  agreement  with  observations.  Second,  the 
longitudinal  correlation  length  of  <^(1)  v^(2)>  Is 


L 


1 

v*<b2> 


[75] 


where  v  is  the  particles  number  density  and  x<b2>  is  the  average  parti¬ 
cle  cross  section.  Finally,  the  transverse  correlation  scale  Is 


Jt  - 


% 


L, 


[76] 


which  is  generally  much  greater  than  the  average  particle  radius  <b> 
characterizing  the  index  covariance  function  <u(l)  u(2)>.  Therefore, 
in  addition  to  the  nonlinearity,  analysis  of  eq.  73  is  further  compli¬ 
cated  by  the  presence  of  the  three  disparate  scales  L,  1,  and  <b>. 

This  is  particularly  troublesome  in  setting  up  finite  difference  meth¬ 
ods  of  solution.  The  results  are  similar  for  <^(1)  vfe(2)>. 

Because  of  the  presence  of  the  derivatives  with  respect  to  z  in 
eq.  73,  the  quantities  <v^.  v^>,  L,  and  JL  are  not  truly  local  proper¬ 
ties  but  they  include  some  form  of  history  of  the  medium  through  which 
the  wave  has  propagated.  In  particular,  there  exists  a  transition 
region  near  the  boundaries  of  the  scattering  medium.  The  size  of  the 
dispersion  domain  varies  with  the  local  conditions  and  is  of  the  order 
of  L  x  A. 


For  the  backscattering  coefficients,  we  need  the  functions  <v^(l) 
v^(2)>  and  <v^(l)  v^(2)>.  These  are  somewhat  easier  to  determine 
because  they  are  not  propagated.  They  define  the  backscattering  prop¬ 
erties  at  the  location  where  the  backscattering  occurs.  Since  the 
scatterers  are  statistically  independent,  the  correlation  does  not 
extend  very  far.  Transversely,  the  functions  <v.(l)  v.(2)>  and  <v.(l) 
v^(2)>  correlate  over  a  region  of  the  same  size  as  the  particles,  i.e. 
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FIGURE  1  -  Diagram  of  primary  reflected  and  refracted  optical  rays  In 
the  plane  of  Incidence  on  a  spherical  particle  of  Index  n 
surrounded  by  a  medium  of  Index  1. 

I  -  <b>.  [77] 

The  contributions  measured  at  points  more  widely  separated  are  uncorre¬ 
lated  because  they  originate  from  different  and  uncoupled  particles. 
Longitudinally,  the  average  thickness  of  the  scattering  sheet  over 
which  there  is  coherence  is  approximated  by  the  mean  distance  between 
the  particles,  l.e. 


1/3 


The  variances  <v^.  and  <v^.  v^>  are  inferred  from  the  propa- 
gated  Vj.  %  and  %b>  by  applying  the  refraction-reflection 

boundary  conditions  at  the  interface  between  the  particles  and  the 
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surrounding  medium.  Figure  1  Illustrates  the  situation  for  the  special 
case  of  spherical  particles*  There  are  two  principal  contributions 
characterized  by  one  external  reflection  and  by  one  refraction  followed 
by  one  internal  reflection,  respectively.  The  incidence  and  exit  an¬ 
gles  defined  in  Fig.  1  are  related  as  follows 


2  9, 

[79] 

2  9  -  2  y. 

[80] 

where 


sin  v  -  -  sin  9  cos  9  [n2  -  2  sin29] 

„2  P 


- 2_  sln  0  ^2  _  gin 2g  [1-2  sin29],  [81  ] 

n2  p 

P 

and  where  n^  is  the  real  part  of  the  particle  refractive  index.  The 
two  main  components  of  <v^.  v^>  (or  <v^.  v^>)  for  a  single  particle  are 
proportional  to  the  variances  <#2>  and  <?2>  calculated  over  the  forward 
scattering  angle  y.  Since  the  scattering  angle  standard  deviation  is 
small,  we  assume  that  y  is  distributed  in  such  a  way  that  the  major 
contributions  to  <$2>  and  <f2>  come  from  small  values  of  y ,  which 
allows  the  linearization  of  eq.  81.  Moreover,  to  model  the  average 
backscattering  effects  of  the  particles,  we  multiply  <$2>  and  <Y2>  by 
the  reflection  coefficient  and  by  the  fractional  cross-sectional  area 
of  the  particles  over  the  scattering  sheet  of  thickness  L  .  He  thus 
obtain 


V, 


%• 


^f> 


(np+i)2 


2/3  9 

v  <b2>. 


[82] 
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<v  v2 

I6n2(2-n  )2  2/3 

a  <V-.  V.  >  - *- - -  V  <b2>, 

<n+l>’ 

[83] 

V  V2 

n2  2/3 

*  <v.  .  V,  >  - 2 -  v  <b2>, 

■*  "*  (n  +1 ) 2 

[84] 

%>■  V. 


16n2(2-n  )2  2(1 

<v,_.  v,  >  - P  ■  P —  v  <b2>. 


~b  ~b 


V1)! 


[85] 


Equations  77,  78  and  82-85  constitute  an  order-of-magnitude  representa¬ 
tion  of  the  covariance  functions  <v^(l)  v^(2)>  and  <v^(l)  v^(2)>. 


Some  of  the  backscattering  coefficients  depend  on  the  correla¬ 
tion  function  H(z,t),£)  defined  by  eq.  72.  Assuming  Gaussian  statistics 
for  the  random  phase  front  fluctuations  ^  and  which  is  a  justified 
standard  approximation,  we  obtain  from  eq.  72 


-  e 


-k2 

Hr  <?f  *  v 


[88] 


where  and  are  the  forward  and  backward  phase  structure  functions. 
For  optical  and  infrared  waves,  the  wave  number  k  is  very  large  and  eq. 
86  can  be  satisfactorily  approximated  by  the  first  term  of  the  series 
expansions  of  P^  and  P^  for  small  separation  distances.  These  can  be 
derived  from  the  solutions  for  <^(1)  vf(2)>.  We  thus  obtain  the  fol¬ 
lowing  transverse  and  longitudinal  correlation  scales  for  H(z,t),£>: 


[87] 


[88] 


a 


t 


■n 
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where  \  is  Che  wavelength.  The  scales  and  L^  are  small  but,  for 
most  practical  applications,  they  remain  much  greater  than  the  respec¬ 
tive  scales  A  and  L  of  <v^(l)  v^(2)>  and  <v^(l)  v^(2)>.  Therefore, 
we  may  set 

H(x.t).fi)  *  1  [89] 

In  eqs.  64-67  with  negligible  consequences  for  the  accuracy  of  the  back- 
scattering  coefficients. 


The  scattering  and  diffusion  coefficients  also  depend  on  the 
amplitude  correlation  functions  6^  and  defined  In  eqs.  70  and  71. 


For  the  coefficients  given  by  eqs.  60,  61,  66  and  67,  the  effects  of  G, 


and  are  negligible  since  the  scale  <b>  of  the  covariance  function 


^r(l)  w(2)>  and  the  scales  A  and  L  of  <v^(l)  v^C2) >  and  <v^(l)  v^(2)> 


are  much  smaller  than  the  scales  of  either  G^  and  G^.  Consequently, 
the  integrals  can  be  calculated  with  good  accuracy  by  substituting  for 
Gj  and  G^  their  values  at  the  origin  which,  by  definition,  are  equal  to 


unity.  However,  in  eqs.  68  and  69,  the  scales  of  and  Gfe  are  compara¬ 


ble  to  those  of  the  phase  front  covariance  functions.  But,  even  in 

these  cases,  the  major  contributions  to  the  integrals  come  from  the 

neighbourhood  of  the  origin  In  the  (t),  p)-space  because  of  the  effect  of 

1  lk  o  ^ 

the  rapidly  oscillating  term  —  exp  [  y  ].  Nevertheless,  the  complete 


tj  *  l2t| 

formal  treatment  of  the  diffusion  coefficients  and  requires  the 


solution  of  G,  and  G. . 

I  D 


Equations  for  G_  and  G,  could  be  derived  from 

£  D 


the  stochastic  eqs.  5-8  but  this  would  produce  multidimensional  differ¬ 
ential  equations  and  more  unknown  moments  for  which  additional  constitu¬ 
tive  relations  would  have  to  be  established.  Instead,  it  is  hoped  that 
the  following  solutions,  derived  by  neglecting  diffraction,  are  suffi¬ 
ciently  accurate  to  account  for  G.  and  G.  In  the  region  adjacent  to  the 


UNCLASSIFIED 

26 


Gf(z,n,£) 


exp 


.  z+tj/2  z+t)/2 
{ - g-  Jdz'/dz"  lia 

O  0  £  -K) 


7  V  : 

~p  ~p 


.  z-q/2  z-rj/2 

-  i  Jdz-  Jdz"  11.  7  7  :  vf(«-.tre/2)> 

O  O  £  +0  K  K 


z+rj/2  z-rj/2 


/dz*  /dz 


lia  7  7  : 

£  *»~P  ~t’ 


<Zt(i'’^^l2)  %s<*' •ve/2’>> 


.  z-rj/2  z+tj/2 

-  i  /dz-  /dz-  11.  V  V  :  ^(z'.r+e/Z)  y^z-.r-g/Z)* 

O  O  £  +0  y  K 


+  i  /dz’  Jdz"  11a  7  7  :  <v^(z'  ,r+£/2)  vf(z",r-£/2)>}.  [90] 

O  O  £  -*0  y  K 


and  a  similar  equation  for  G^(z,t)»£)>  This  will  have  to  be  verified 
experlaentally  but  the  errors  are  expected  to  be  small  since  the  prin¬ 
cipal  contribution  to  the  Integrals  of  eqs.  68  and  69  which  coae  froa 
the  neighbourhood  of  the  origin  are  independent  of  the  approximation. 
The  alternative  is  a  auch  higher  level  of  mathematical  coaplexity  or 
the  need  for  some  empirical  input. 


Equations  62-69  and  89-90  give  the  propagation  coefficients  in 
terms  of  the  forward  and  backward  phase  front  covariance  functions. 
These  functions  satisfy  nonlinear  partial  differential  equations  which 
were  not  solved  in  this  report.  The  principal  difficulties  are  the 
nonlinearity  of  the  equations  and  the  presence  of  three  dlsparste 
scales:  the  average  particle  size,  the  aean  distance  between  wave 
interactions  with  the  particles,  and  the  latter  multiplied  by  the 
forward  scattering  phase  front  angle  standard  deviation.  What  was 
achieved  in  this  report  is  the  deteralnatlon  of  these  scales  and  the 
derivation  of  order-of-aagnltude  estimates  of  the  pertinent  variances. 
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Consequently,  the  present  state  of  the  node I  allows  only  rough  esti¬ 
mates  for  the  propagation  coefficients. 

Although  a  satisfactory  solution  of  the  nonlinear  eq.  73  has  yet 
to  be  derived,  eqs.  58-69  together  with  eq.  73  for  the  phase  front 
angle  covariance  function  provide  a  formal  framework  for  relating  the 
propagation  coefficients  to  the  medium  properties.  Moreover,  if  eq.  90 
holds,  this  problem  is  totally  uncoupled  from  the  propagation  problem 
since  the  wave  amplitude  is  not  in  eqs.  58-69.  Therefore,  the  determi¬ 
nation  of  the  scattering  and  diffusion  coefficients  is  a  well-posed 
problem,  albeit  difficult,  that  can  be  handled  independently  of  propa¬ 
gation.  The  situation  is  somewhat  similar  to  the  more  common  approach 
where  Mie  calculations,  needed  to  specify  the  phase  function,  are 
performed  independently  of  the  solution  of  the  radiative  transfer 
equation. 

Equations  58-69  are  valid  for  inhomogeneous  media,  the  only 
restriction  being  that  the  size  of  the  lnhomogeneltles  must  be  larger 
than  the  beam  diameter.  It  must  also  be  noted  that  the  various  scat¬ 
tering  and  diffusion  coefficients  are  not  truly  local  values.  They 
bear  some  history  of  the  scattering  medium  over  a  region  defined  by  the 
convolution  of  the  coherence  domain  of  the  phase  front  angle  covariance 
functions  with  the  integral  operators  0*  and  0**. 

The  physical  interpretations  are  straightforward.  The  scat¬ 
tering  coefficients  are  responsible  for  the  transfer  of  energy  from  the 
coherent  state  into  the  incoherent  or  scattered  state,  and  from  the 
forward  propagating  beam  into  the  backward  propagating  beam  and  vice 
versa.  The  diffusion  coefficients  measure  the  rate  at  which  the  scat¬ 
tered  energy  is  spatially  diffused,  or  the  rate  at  which  the  scattered 
beam  spreads  out.  They  account  for  the  multiple  scatterings. 

The  diffusion  coefficients  are  dyadics  or  second-order  tensors 
but  all  the  components  are  not  independent.  From  Hypothesis  5,  the 
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phase  front  angle  covariance  functions  are  statistically  homogeneous 
and  Isotropic  In  the  plane  normal  to  the  propagation  axis*  The  general 
form  for  such  an  Isotropic  tensor  is  (Ref.  21) 

<v(l)v(2)>  -  M(z,  t),  p2)  £g  +  N(z,T),p2)  £,  [91  ] 

where  6  is  the  two-dimensional  unit  dyad  or  the  two-dimensional 
Kronecker  delta  and  M  and  N  are  scalar  functions.  The  diffusion  coef¬ 
ficients  are  thus  given  by  an  expression  of  the  form 

iko2 

z  &  ( Z  ® 

-  /  -1  e  f  //  d2fi  e  11  [M(z,h»P2)je£  +  N(*,t|,p2)  6] 

z  •  lkp2 

"  "  TT  \  T*  «‘Af(Z,Tl)/  Pdp  ®  2r>  [M(z,twp2)+2N(z,ti,p2)]  6,  [92] 

or  in  shorter  notation  by 

»-»§.  [w] 

where  D  is  a  scalar  function  of  position  z  along  the  direction  of 
propagation. 

The  mathematical  difficulty  of  solving  exactly  for  the  scat¬ 
tering  and  diffusion  coefficients  will  be  dealt  with  in  a  subsequent 
work.  In  this  report,  we  will  proceed  with  the  solution  of  the  propa¬ 
gation  equations.  Comparison  with  experimental  data  or  empirical 
modeling  can  always  be  used  to  determine  the  coefficients. 


The  governing  propagation  equations  are  obtained  by  substituting 
the  constitutive  relations,  eqs.  40-52,  for  the  higher  order  unknown 
moments  In  eqs.  23,  24,  27  and  28.  We  limit  the  present  analysis  to 
the  more  practical  case  where  there  Is  only  one  coherent  beam  transmit¬ 
ted  In  what  Is  chosen  as  the  forward  direction.  Under  this  condition 
<A^>  ”  0.  Moreover,  we  assume  that  the  aerosol  diffusion  effect  Is 
much  more  Important  than  the  diffraction  spreading  of  the  beam,  which 
is  well  justified  In  most  applications.  We  can  thus  neglect  the  terms 
of  the  form  .  [<a*7^a>  -  <aVja*>].  With  these  approximations  and 

using  eqs.  36  and  37,  we  obtain  the  following  set  of  propagation 
equations: 

a)  for  the  average  reduced  amplitude  A  ■  <A  >: 
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aa  -  2k<W>.  [97] 

Equations  94-96  constitute  the  proposed  propagation  model*  They 
form  a  closed  set  of  three-dimensional,  coupled,  partial  differential 
equations  for  the  solution  of  the  forward-  and  backscattered  lrradlance 
beam  profiles  resulting  from  a  coherent  optical  or  infrared  beam  di¬ 
rected  into  a  particulate  medium.  The  dimensions  are  reduced  to  two 
under  cylindrical  symmetry. 

The  various  absorption,  scattering  and  diffusion  coefficients 
are  formally  related  to  the  medium  properties  through  eqs.  58-69. 
Homogeneity  is  assumed  only  on  the  scale  of  the  beam  diameter.  The 
model  handles  all  orders  of  scattering  in  a  global  fashion  but  cannot 
differentiate  between  them.  However,  the  single-scattering  results  can 
be  derived  by  dropping  eqs.  95  and  96  and  solving  eq.  94.  Similarly, 
the  first-order  multiple-scattering  approximation  for  the  backscattered 
beam  can  be  modeled  by  setting  I*  ■  0  in  the  right-hand  side  of  eq.  96, 
and  by  arbitrarily  adding  the  forward-scattering  coefficient  Re(acf)  to 
[ofl  -  Re(0a|,8)  +  Re(<rfe)]  *n  t*ie  left-hand  side  of  the  same  equation  to 
account  for  the  single-scattering  losses. 

It  is  emphasized  that  the  model  predicts  the  lrradlance  and  not 
the  radiant  Intensity.  Experimentally,  this  means  that  (A^AJ  +  I+)  and 
I~  must  be  measured  with  wide  field  of  view  detectors. 

6.0  BEAM  WAVES 


We  consider  the  problem  of  an  initially  Gaussian  beam  illuminat¬ 
ing  a  particulate  medium.  The  propagation  coefficients  are  either 
specified  directly  or  have  been  derived  from  eqs.  58-69.  Then  the 
propagation  eqs.  94-96  must  be  solved.  The  principal  difficulty 
resides  in  the  coupling  between  the  forward-  and  the  backscattered 
lrradlance.  In  the  present  application,  the  coupling  is  treated 


UNCLASSIFIED 

31 


iteratively,  first  by  setting  I"  •  0  in  eq.  95  and  solving  for  1+  and 
I~  in  sequence,  and  then  by  substituting  the  I~  solution  back  into  eq. 
95  and  repeating  the  process.  Under  these  conditions,  the  equations 
can  be  solved  analytically. 


To  reduce  the  number  of  parameters,  we  assume  a  collimated  beam, 
i.e.  F-*»,  and  we  make  the  following  simplifications 


a  , 
af 

*  a  , 
ab 

“  °af  s 

“  aabs' 

<r. 

a  cf. 

a  (T 

»  (f 

sf 

sb 

cf 

cb  * 

[98] 

[99] 


which,  from  the  defining  eqs.  58-67,  are  expected  to  hold  for  particles 
larger  than  the  wavelength.  Furthermore,  we  define  the  nondimens ional 
independent  variables 


z 


T)  =  /  ot(2,)dz’, 
o 

[100] 

P  =  r/wo. 

[101] 

where 


a  .  +  a  ,  +  cT, 
af  cf  cf 


]. 


[102] 


and  where  w  is  the  initial  beam  radius  at  1/e  in  irradiance.  The 
o 

nondimens ional  propagation  equations  thus  become 


OA 

^T  +  T<1  +  !*>*,-“  7pAr-°- 

+  a  1+  -  K.  V2  1+  -  0  A  A*  +  x  I“, 
or)  f  p  r  r  r  ’ 


[103] 


[104] 
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+  a  1“  -  <w  V2  r  -  x  A  A*  +  x  I4 . 
dti  bp  r  r 

The  similarity  parameters  are 


[105] 


I«  [v  cs[  +  act  +  <£f] 


A  -  1/(2  t»2  <,t). 


"tV  °afs+  Vi 
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[106] 
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[108] 


„  RetV+  “cfl 
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Re[Df] 
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°tWo 


[109] 

[no] 

[in] 

[112] 


For  the  calculation  results  presented  here,  we  have  taken  x»  a, 
p  and  x  to  be  constant  throughout  the  medium.  This  constitutes  a  rea¬ 
sonable  approximation  if  only  the  particle  number  density  is  allowed  to 
change  with  little  variations  in  the  nature  of  the  particles  and  their 
size  distribution.  The  parameter  A  varies  with  a but  it  is  generally 
too  small  to  make  appreciable  differences  in  the  solutions  and  it  was 

assumed  constant.  Moreover,  the  coefficients  k r  and  k.  increase  from 

r  d 

zero  to  asymptotic  values  in  their  respective  transition  regions  near 
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the  medium  boundaries.  However,  since  the  scattered  irradiance  is 
saall  in  the  transition  regions,  we  have  approximated 

*  icb  -  tc  ■  cat  [ll3] 

everywhere  with  little  effect  on  the  resulting  solutions. 

The  solutions  were  calculated  for  the  boundary  conditions 


Ar(o,p)  -  exp  (-p2),  [114] 

I+(o,p)  -  0,  [115] 

I'(C.P)  “  0,  [116] 

where  C  is  the  optical  depth  of  the  cloud,  l.e. 

Z 

C  -  /  ot  (z’)dz'.  [117] 

o 


The  solutions  were  calculated  analytically  and  only  the  first  Iteration 
was  worked  out.  The  results  are  presented  in  Figs.  2-5  for  the 
following  set  of  parameters:  p  ■  0,  A  ■  4.7  x  10-5,  a  *  i  *  0.03,  p  ■ 
0.97  and  k  ■  65.  These  are  the  estimated  similarity  parameters  for  a 
1.06-  |in,  0. 15-m-diameter  laser  beam  propagating  in  an  advective  fog  of 
mass  loading  equal  to  0.5  g(H20)/m3  with  a  modal  particle  radius  of 
5  pa.  The  water  droplets  were  assumed  to  be  nonabsorbing. 

The  transmitted  beam  irradiance  profile,  i.e.  (A^A*  +  I+),  at 
T)  ■  0,2, 4, 6  and  8  optical  depths  are  plotted  in  Fig.  2.  The  irradiance 
is  normalized  to  unity  on  the  propagation  axis  (p  ■  0)  at  the  transmit¬ 
ter  (tj*  0).  The  results  show  that  the  basic  Gaussian  shape  of  the 
beam  is  not  appreciably  broadened  up  to  optical  depths  as  large  as  6. 
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FIGURE  4  -  Predicted  Integrated  backacattered  lrradiance  from  an  originally  Gaussian  beam 
propagated  through  a  cloud  of  optical  thickness  ;  as  a  function  of  c*  Simi¬ 
larity  parameters  as  In  Fig.  2.  Upper  curve  is  result  of  present  multiple¬ 
scattering  model  while  the  lower  curve  represents  predictions  from  first-order 
multiple-scattering  version  of  the  model. 
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Whereas  the  effects  of  nultiple  scattering  were  unimportant  for 
on-axis  measurements  of  the  transmitted  beam.  Fig.  4  shows  that  they 
are  significant  for  the  analysis  of  the  detected  backscatter.  This  is 
even  more  dramatic  if  we  consider  the  lidar  signal,  i.e.  the  range¬ 
gated  backscatter  signal,  as  illustrated  in  Fig.  5.  At  optical  depth 
3,  the  difference  between  the  full  multiple  scattering  and  the  first- 
order  multiple-scattering  predictions  has  already  reached  a  multiplica¬ 
tive  factor  of  10  and  it  continues  to  grow  with  increasing  optical 
depth.  These  results  have  great  implications  in  the  interpretation  of 
lidar  measurements.  Careful  comparison  with  well-documented  data  is 
still  needed  but  Figs.  4  and  5  indicate  that  the  effect  of  multiple 
scattering  is  potentially  very  important. 

The  calculation  results  depicted  in  Figs.  2-4  agree  well  with 
various  known  effects  of  the  multiple-scattering  phenomenon,  in  parti¬ 
cular  its  greater  influence  on  the  backscatter  than  on  the  transmission 
measurements.  These  solutions  were  derived  analytically  with  the  help 
of  a  number  of  simplifying  approximations  stated  at  the  beginning  of 
this  chapter.  For  more  general  applications,  the  propagation  eqs. 
103-105  can  be  solved  numerically  with  no  restrictions  on  the  absorp¬ 
tion,  scattering  and  diffusion  coefficients  other  than  the  hypothesis 
of  homogeneity  over  a  scale  of  the  order  of  the  beam  diameter.  The 
equations  present  no  particular  difficulties  for  a  numerical  solution. 

7.0  PLANE  WAVES 


It  is  Instructive  to  consider  the  special  geometry  of  plane  or 
spherical  waves.  This  simplification  allows  complete  solutions  to  be 
derived  in  closed  analytic  form.  These  analytic  expressions  are  very 
useful  for  parametric  studies.  Moreover,  the  planar  or  spherical  geom¬ 
etry  represents  adequately  many  physical  situations,  e.g.  the  planetary 
atmospheres,  the  clouds  and  the  ocean  waters  illuminated  by  the  sun. 
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The  sene  simplifying  approximations  as  In  the  preceding  section 
are  made*  The  nondimenslonal  propagation  eqs*  112-114  are  adapted  to 
planar  geometry  by  simply  neglecting  the  transverse  gradients.  We  thus 
obtain 


dA  A* 

-/  r  +  A  A*  -  0, 
dii  r  r  * 


Hf  +  «  ^  "  Mr A*  +  x  I", 


[118] 

[119] 


dl  ^  T- 
•x—  +  a  I 
dn 


x  A  A*  +  x  I+. 
r  r 


[120] 


Equations  118-120  form  a  closed  set  of  three  simultaneous,  coupled 
ordinary  differential  equations.  The  boundary  conditions  are 


ArA*  (0)  -  1, 

[121] 

1+  (0)  -  0, 

[122] 

I"  (O  -  0, 

[123] 

where  C  Is  the  optical  thickness  of  the  plane-parallel  scattering 
medium. 


The  set  of  eqs.  118-120  is  identical  to  the  four-flux  theory 
described  In  Chapt.  10  of  Ref.  1  where  A^A*  is  the  forward  collimated 
flux  and  1+  and  I“,  the  forward-  and  backward-diffused  fluxes  respec¬ 
tively.  We  have  simplified  the  problem  by  assuming  no  backward- 
collimated  flux  and  making  such  approximations  as  defined  by  eqs.  98 
and  99.  These  make  the  algebra  easier  but  they  are  not  fundamental  to 
the  model.  Therefore,  the  four-flux  theory  is  the  plane-wave  asymp¬ 
totic  form  of  the  propagation  model  proposed  in  this  report.  However, 
contrary  to  the  four-flux  theory,  which  is  constructed  heuristlcally 
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with  propagation  coefficients  not  clearly  related  to  the  physical  pa- 
raaeters  of  the  particulate  aediua,  our  model  Is  derived  froa  Maxwell's 
equations  and  the  coefficients  are  foraally  related  to  the  aediua 
properties. 


The  propagation  eqs.  118-120  with  the  boundary  conditions  121- 
123  are  straightforward  to  solve.  The  siailarlty  parameters  are  the 


M°c£  +  »rf] 


[124] 


and  the  ratio  of  the  backscattering  to  the  total  scattering 


coefficients 


[125] 


With  the  further  approxiaatlon  ogf  «  acf,  which  is  almost  always  sat¬ 
isfied  in  practical  cases,  the  results  are: 

a)  for  the  total  lrradlance  transmitted  through  a  cloud  of 
optical  thickness  C 


A  A*(C)  +  I+(C)  -  - — - x-r  , 

(6fe)  -  (6-e)e-ze‘’ 


[126] 


and  b)  for  the  integrated  backscattered  lrradlance  froa  a  cloud  of 
thickness  C 


r<o,o  - 


no  [i 


-  .-2cC 
6 


(W*e)  -  ( 6— e>  e 


-2eC  * 


[127] 
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where 

6  -  (1-fl)  +  DO,  [128] 

e*»  (1-0)  [(1-0)  +  2IIQ]. 

The  corresponding  single-scattering  (SSCT)  and  first-order  multiple- 
scattering  (FOMS)  solutions  are 


[130] 

t1'*0-'*  W  -  ®  [1 '  **2e)- 

[131] 

For  spherical  waves,  we  find  exactly  the  same  expressions  except 
for  the  multiplicative  geometric  factors  R22  and  R2^  on  the  left-hand 
side  of  eqs.  126(130)  and  127(131)  respectively.  Rj  and  R2  are  the 
radial  positions  of  the  centered  front  and  back  spherical  boundaries  of 
the  scattering  medium. 

The  transmitted  lrradiance  solutions  are  plotted  in  Fig.  6 
versus  the  cloud  optical  thickness  C»  for  II  ■  0.04  and  various  values 
of  the  albedo  Q.  The  decay  rate  with  respect  to  C  varies  inversely 
with  0  but,  in  all  cases,  it  remains  smaller  than  for  the  single- 
scattering  solution  (SSCT).  The  latter  is  reached  in  the  limit  0*0, 
l.e.  when  the  absorption  becomes  dominant  over  the  scattering. 

The  integrated  backscatter  results  are  displayed  in  Fig.  7.  All 
the  curves  saturate  to  some  asymptotic  level  but  the  required  cloud 
thickness  and  the  saturation  level  increase  with  Q.  If  there  is  no 
absorption  (Q  ■  1),  the  transmitted  lrradiance  is  eventually  all 
reflected  back  for  optical  thicknesses  much  greater  than  X/n.  In  the 
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FIGURE  7  -  Predicted  Integrated  backscattered  lrradlance  In  a  plane-parallel  geometry. 

Notation  aa  In  Fig.  6  except  for  curve  marked  FOMS  Which  represents  the  result 
of  flrat-order  multiple-scattering  approximation. 
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first-order  multiple-scattering  limit,  the  backscattered  irradiance 
saturates  quickly  for  C  *  2  and  at  a  level  equal  to  IIQ/2  which  repre¬ 
sents  only  a  small  fraction  of  the  transmitted  signal. 

Figures  6  and  7  show  large  differences  between  the  plane-wave 
multiple-  and  single-scattering  solutions.  This  is  in  agreement  with 
the  fundamental  differences  between  the  two  approaches.  In  single 
scattering,  the  scattered  energy  is  assumed  lost  whereas  it  is  simply 
redistributed  spatially  in  multiple  scattering.  Because  of  the  plane- 
parallel  geometry,  the  losses  at  one  point  due  to  spatial  redistribu¬ 
tion  are  compensated  for  by  gains  from  the  adjacent  regions.  Hence, 
the  multiple-scattering  solutions  should  indeed  predict  higher  trans¬ 
mitted  and  backscattered  irradiance  levels,  as  illustrated  by  our 
results  in  Figs.  6  and  7.  Moreover,  it  is  well  understood  that  the 
single-scattering  solution  should  become  more  accurate  as  the  effect  of 
absorption  increases  since  the  absorption  represents  a  true  energy  loss 
to  propagation.  This  is  also  well  verified  by  the  predictions  of  our 
model . 


An  interesting  physical  situation  is  the  no-absorption  limit, 
i.e.  Q>1.  In  this  case,  6-*H,  e-*0,  and  the  transmitted  and  integrated 
backscattered  lrradiances  become  respectively 


ArA*(C)  +  I+(C>  "  i+Hc  '  tl32j 

^(O.O-T^C*  U33l 

Bucher  (Ref.  16)  has  shown  that  Monte  Carlo  simulations  fit  very  well 
the  expression 


It<0 


1.69 _ 

1.42  +  (l-<cose>)C  ’ 


[134] 
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for  (1-<cos0>)C  >  1*  where  I£  ■  AfAj  +  1+  and  where  <cos  0>  is  the 
average  cosine  of  the  scattering  angle  over  the  complete  forward  and 
backward  directions  relative  to  the  incident  axis.  It  is  not  clear  how 
(l-<cos  0>)  relates  exactly  to  our  parameter  n,  but  they  both  vary  in 
the  same  direction  with  the  ratio  of  the  backseat ter ing  to  the  forward¬ 
scattering  coefficients  and  they  are  both  of  the  same  order  of  magni¬ 
tude,  i.e.  a  few  percent.  Therefore,  eqs.  132  and  134  show  that  our 
analytic  solution  is  in  good  qualitative  agreement  with  Monte  Carlo 
data.  This  functional  relationship  was  also  verified  in  a  laboratory 
experiment  (Ref.  23). 


8.0  CONCLUSION 


We  have  derived  from  Maxwell's  equations  a  propagation  model  for 
optical  and  infrared  waves  in  particulate  media.  It  has  the  form  of 
simultaneous  three-dimensional  partial  differential  equations  for  the 
average  amplitude  and  the  average  scattered  irradiance  of  the  forward- 
and  backward-propagating  waves.  These  equations  involve  absorption, 
scattering  and  diffusion  coefficients  that  are  formally  related  to  the 
statistical  properties  of  the  particles  and  the  surrounding  medium. 

The  model  handles  beams,  includes  all  orders  of  scattering,  and  is 
applicable  to  inhomogeneous  media  provided  the  inhomogeneity  scale  is 
not  smaller  than  the  beam  diameter.  The  main  restriction  is  that  the 
fraction  of  volume  occupied  by  the  particles  to  the  $th  power  must  be 
much  smaller  than  unity,  which  is  well  verified  in  atmospheric  applica¬ 
tions.  The  polarization  effects  are  neglected  in  the  present  version. 

The  propagation  equations  are  linear  and  simple  to  solve,  at 
least  numerically.  The  coefficients  are  related  to  the  medium  proper¬ 
ties  through  nonlinear  equations  which  are  difficult  to  deal  with.  In 
this  report,  only  order-of -magnitude  estimates  of  these  coefficients 
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were  obtained.  The  important  observations,  however,  are  that  the  rela¬ 
tions  do  exist  and  that  they  are  not  coupled  to  the  propagation  prob¬ 
lem.  Thus,  the  coefficients  can  be  solved  independently  or  inferred 
by  some  other  means  such  as  experimental  measurements  or  empirical 
modeling. 

The  solutions  derived  in  this  report  are  in  good  agreement  with 
the  present  understanding  of  the  multiple-scattering  phenomenon. 

Closed  form  expressions  were  obtained  that  are  useful  for  analytical 
purposes.  In  particular,  the  model  provides  a  convenient  theoretical 
tool  for  the  analysis  of  the  multiple-scattering  effect  on  the  inver¬ 
sion  of  the  lidar  data.  This  will  constitute  the  subject  of  a  future 
report. 
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